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Abstract. A XY Heisenberg spin chain model with two perpendicular spins par 
site is mapped onto a Kirchhoff thin elastic rod. It is shown that in the case of 
constant curvature the Euler-Lagrange equation leads to the static sine-Gordon 
equation. The kink-antikink type and periodical static solutions for these models 
are derived. 
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1. Introduction 

The study of an elastic rods pQ is a subject to increased interest [3 El EI El El El EJ 
especially in connection with the biomathematical models of proteins and of DNA 
[HI E31- The main feature of a thin rod is a space curve (rod's axis) and the 
corresponding orthonormal frame with a tangent vector t to the axial curve. The 
static energy of the elastic rod is related to the bending and twisting energies. It is 
tempting to map the elastic rod problem to a classical spin chain 11 (in the continuum 
limit, where the normalized spin S is mapped onto the tangent vector t). We will show 
however that the full mapping of the elastic rod onto a spin-chain model requires a 
system of two orthogonal spins. 

The spin Hamiltonian for a Heisenberg spin chain is given by the following 
expression: 

H = J J2s i (s)-S i+1 (s), S?=S? +1 = 1. (1) 

i 

In the continuum limit this Hamiltonian goes over to 

f +oc fdS(x)\ 2 

In the case of XY spin chain the spin is given by the rotation angle 6(x): S(s) = 
(cos 6(x), sin 0(x)). The Hamiltonian now reads H = J J^°° (|§) 2 dx. 
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This letter is organized as follows: In section 2 are discussed the static properties 
of the Kirchhof equations for a thin elastic rod. It is shown that if a curvature is 
present the twist angle satisfies the static sine-Gordon equation. The mapping of the 
Kirchhoff model onto a XY spin chain is done in Section 3. In Section 4 the soliton-like 
solutions for the static sine-Gordon equations are briefly discussed. 



2. Kirchhoff Model for Elastic Rods 



k(s) 



The model introduced by Kirchhoff (1859) describes the shape and the dynamics 
of a thin elastic rod in equilibrium and is based on the analogy with the dynamics 
of a heavy spinning top (the Lagrange case). The shape is described by the static 
Kirchhoff model while the time evolution - by the dynamical Kirchhoff model. Here 
we will concentrate on the statics of thin elastic rods (here and below we shall call 
them Kirchhoff rods). 

We consider a static space curve R(s) : R — > M 3 as a smooth function mapping 
the arc-length interval J Cl into the physical space IR 3 . For every s we define the 
Frenet basis (t(s), n(s), b(s)) to be the normal, binormal and the tangent vectors to 
the curve(s). The tangent vector is a unit vector given by t = (^p) and the curvature 
k(s) of the curve at the point s is then given by: 

dt 

ds 

The triad (t(s), n(s), b(s)) evolves along s according to the Frenet-Serret equations: 

dt dn . . . . dh 

- = .n( s ) - = - K t( s ) + rb(s) ^ = -m(s), (3) 

where r is the torsion of the curve R(s). If the curvature k and the torsion r are 
known for all s then the Frenet-Serret triad can be obtained as unique solution of 
Next the space curve R(s) can be reconstructed by integrating the tangent vector 

A thin rod can be modelled by a space curve R(s) joining the loci of the centroids 
of the cross sections together with the local basis (di(s), d 2 (s), d 3 (s)) attached to the 
rod material. This local basis can be expressed through the Frenet-Seret triad as 
follows: 

/ I 

(d 3 (s),d 2 (s),di(s)) = (t(s),n(s),b(s)) cos</> -sin<; 

y sin0 coscj) 

where 4> is the twist angle of the rod. The components of the derivatives of the local 
basis (d 3 (s), d 2 (s), di(s)) with respect to s can be expressed by using the twist vector 
k(s) = Kidi + K 2 d 2 + K 3 d 3 as follows: 

T=k(s)xd,(s), i = 1,2,3. 
ds 

The static Kirchhoff equations describe the shape of the rod under the effects of 
internal elastic stresses and boundary constraints, in the absence of external force 
fields. Let F(s) is the tension and M(s) is the torque of the rod. In the approximation 
of a linear theory (the Hook's law applies) the torque M is related to the twist vector k 
by M(s) = S ■ k(s), where S — diag (1, a, b). The constant a measures the asymmetry 
of the cross section and b is the scaled torsional stiffness. In particular for symmetric 
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(a = 1) hyperelastic (b = 1) rods we have M(s) = k(s). In the generic case the torque 
is 

M(s) = K 1 (s)d 1 (s) + an 2 (s)d 2 (s) + bK 3 (s)d 3 (s) (4) 
and the elastic energy of the Kirchhoff rod is given by: 

H= \ f 2 M(s) ■k(s)ds = \ f 2 (kI{s) + aK 2 2 (s) + bK 2 3 (s)) ds (5) 

^ J si ^ J si 

The conservation of the linear and angular momenta is provided by the static Kirchhoff 
equations: 

£-0, (•) 
ds 

? + d 3 ( S )xF( s ) = 0. (7) 
ds 

Here F(s) is the tension of the rod and the torque M(s) is given by Q and the twist 
vector reads 

k(s) = (k(s) sin <j), K (s) cos (j), t + (j) s ), (8) 

where the twist angle is a function of the arc length parameter s: ef> = (f>(s). The 
expression for the tension F(s) in the local basis di, da, d3: 

F(s) = F 1 (s)d 1 (s) + F 2 (s)d 2 (s) + F 3 (s)d 3 {s) 
reduces the Kirchhoff equations © and to the following system of ODE's: 

F hs + K2F 3 ~ k 3 F 2 = (9) 

F 2 . s + k 3 ^i - kiF 3 = (10) 

^3,, + k,F 2 ~ k 2 F x = (11) 

F x = -ok 2 , s + (b - 1)«ik 3 (12) 

^2 = «i,s + (b - a)n 2 K 3 (13) 

bn 3:S + (a - 1)kxk 2 = 0. (14) 

Using the parameterization of the twist vector © from (|14|l for the case of constant 
curvature k(s) = kq we get the famous static (scalar) sine-Gordon equation: 

(a_l) 
ds 2 ' b 

This second order differential equation is a completely integrable Hamiltonian system 
and allows so-called "soliton"-like solutions. It appears in a wide variety of physical 
problems for e.g. charge-density-wave materials, splay waves in membranes, magnetic 
flux in Josephson lines, torsion coupled pcndula, propagation of crystal dislocations, 
Bloch wall motion in magnetic crystals, two-dimensional elementary particle models 
in the quantum field theory, etc. 

For a symmetric rod, i.e. a = 1 the equation l|15|l simplifies to = 0, i.e. 
cf) s = const. This is the usual case widely discussed in the literature [T5]. The cross- 
section of the symmetric rod (a = 1) has a continuous rotational symmetry around 
the central axis. Therefore the elastic energy density h does not depend on <j) and from 
the variational principle it follows that h could be only a function of the derivatives 
<p s . Therefore the only solution for a constant twist is <j) s = const. 

There have been even attempts to generalize this result to the asymmetric case 
(o^l) i.e. to show that <j) s = const holds true for any Kirchhoff rod. 



— + u k 2 q smu(s) = 0, u(s) = 2(j>(s). (15) 
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The asymmetric case has been overlooked for a long time. For a constant 
curvature and torsion along the centre-line it represents another integrable case of 
the Kirchhoff equations for a thin elastic rod. This opens new possibilities for a 
more adequate modelling of bio-polymers and gives the phenomenological bases for 
the widely used DNA models. Here there is no more continuous rotational symmetry 
of the cross-section around the centre axis and obviously h depends on (ft as well. So 
the constant twist is no more a solution. 

The solution of 1)15(1 is compatible with the full system of Kirchhoff equations 
© for constant curvature and torsion. 



3. The Spin Chain Model 



Let us consider the following spin chain model with two perpendicular spins per site: 
the spin vectors Si and S2 have different lengths and are given by 

S? = l, <ft = 4>(s); (16) 
S 2 (s)=-\l - " ' * n(s)sin<ft + J b " + ] b(s) cos(ft, 



Si (s) = n(s) cos (ft -I- b(s) sin (ft, 



1 



1 



a 2 — 



1 



Here also both spin vectors are orthogonal:Si(s) • 82(5) = for every s. This is an 
integrable system with Hamiltonian given by 

H = Jo^2 (Si,i(s) ■ S M+ i(s) + S2,i(s) ■ S 2 ,i+i(s)) , 



which in the continuum limit leads to: 



H = Jo 



ggi 

ds 



dS 2 

ds 



ds = J\ 



(6(t + <ft s ) 2 + K 2 (a- l)sin 2 



ds, 



where J\ is the renormalized coupling constant, s s [s±, S2] and the subscript s means 
a derivative with respect to s. This Hamiltonian coincides with that one in JSJ) where 
k(s) has been replaced from JSJ. Thus the asymmetric Kirchhoff rod is mapped onto 
a two-spin chain. If the curvature is constant k(s) = kq, then the Euler-Lagrange 
equation gives the (scalar) static sine-Gordon equation: 

G? 2 <ft o 1 — a 

~rr + K n — ; — sin (ft cos (ft = 0. (17) 
ds z 

When a 1 the Hamiltonian ((17(1 simplifies to (r + <ft s ) 2 ds, so the Euler-Lagrange 




;(D 

i+1 




i+1 



Figure 1. A two-spin XY chain system with a coupling constant Jo that 
corresponds to the Kirchhoff rod model. 
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equation leads to cf) ss — 0, or <j> s = const. Note that the case a = 1 corresponds to 
a symmetrical Kirchhoff rod. Such a rod-model is mapped onto a symmetrical spin 
chain system (S 1; S2) with = S2 = 1. Note that the mapping of the Kirchhoff 
symmetric rod needs a two-spin XY chain rather than a simple one-spin XY chain. 
From 117f) one can easily get that 
2 



,1 



• sin <f){s) = 0, 



ckp 
ds 

which may be used for the calculation of the corresponding static energy. For the 
asymmetric thin elastic rod the model Hamiltonian takes the form: 



H = 



[t + 4> s ) 2 + 4- 



1 



ds 



(18) 



Here we shall discuss in brief the static soliton solutions of the model (|18fl : The kink 
type solution of the static sine-Gordon equation (JTSJ, l(T7|) is given by 



2 arctan exp 
and the corresponding static energy is 
-Ekink = 4kq 




1 -a 



tanh kq 



1 -a. 



(19) 



(20) 



(21) 



with the periodicity 4^r \J K (k) , where k is the modulus of the Jacobian elliptic 

function sn (sine amplitude), and K(k) is the complete elliptic integral of the first 
kind. In the limit k — > 1 we have K(k) — > 00 and the half-period tends to infinity as 
well, so we recover the single kink soliton solution (|19|l . 

The corresponding static energy per soliton of the soliton lattice is given by: 



b I V & 

The periodic (soliton lattice) solution of l|15fl . (|17fl is 

6(s) = 2 arccos 




k 



E(k) - -{k'YK{k) 



(22) 



where E(k) is the complete elliptic integral of second kind. In the single soliton limit 
(k — > 1) the lattice energy per soliton reduces to eqn. lj2U|) . 



4. Conclusions 



We have shown that the single asymmetric elastic Kirchhoff rod model can be mapped 
onto a 2-spin XY Heisenberg chain and the spin vectors must have different lengths. 
In this case the Euler-Lagrange equation for the spin chain Hamiltonian gives the 
static sine-Gordon equation. For the case of symmetric rods (a — 1) both spins have 
the same length. The symmetric (a = 1) and the asymmetric (a / 1) Kirchhoff rods 
have very different static properties. In general the family of thin elastic Kirchhoff 
rods falls into two groups: i) the group of symmetric rods (a = 1). Here the twist is 
constant along the rod and if the torsion is constant as well, the Kirchhoff equations 
are integrable and the curvature satisfies the non-linear Schrodinger equation; ii) the 
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group of asymmetric rods (a ^ 1). Here in general the twist and the curvature satisfy 
a coupled differential equations (for a constant torsion) . In the special case where the 
curvature is constant the system of Kirchhoff equations is integrable again and the 
twist satisfies the sine-Gordon equation. 

The dynamics of such models, which is of interest for realistic biopolymers, should 
be investigated. Due to the Galilean invariance of the sine-Gordon equation 115(1 a 
special class of dynamical travelling wave type solutions can be obtained from the 
static ones by Galilean boost. 

The general assumption that all thin rods exhibit constant twist should now be 
restricted to the class of symmetric thin rods only and to all straight rods as well. 
The class of asymmetric thin rods does not belong to this category. Here the twist is 
not constant and "interacts" with the curvature. In the case of constant curvature, 
the problem has exact solution. For non-constant curvature the case is more complex 
and should be of considerable interest e.g. for the problem of DNA supercoiling 

[H m mm HUGH]. 
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